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1 Introduction 



In ^T] Gould and Hopper generalized the classical Hermite polynomials by 
introducing the polynomials i/^(x,a,p) defined by 

Hn{x,a,p) = {-l^x "exp(^5x'^) — [x" exp (-px*^)] . 

They derived several properties of H^{x,a,p), including a generating func- 
tion, differentiation, addition and operational formulas. As they remarked, 
the case a = was considered before them by Bell in P, where he analyzed 
the polynomials 

{x, t; r) = exp {-xf) — exp {xf ) . 

Dhawan ^ , obtained a new generating function, summation formulas, a 
hypergeometric representation and some integrals of H^{x, a,p), including 

oo 

y^^exp {—px'^) H^{x, a,p)H^{x, a,p)dx = 



for n > ms. 

Kalinowski and Seweryhski jT3], |2I], constructed a differential equation 
of order r for H^{x,0,p). They also proved that the polynomials H^{x,0,p) 
are not orthogonal with respect to the weight function exp {—px^) , r even, 
in the interval (—00,00) . 

Todorov 25J, considered the polynomials H^{x,0,p) in connection with 
his analysis of the n*'' derivative of the composite function / (x^) . Unaware 
of Kalinowski and Seweryhski's previous work, he re-derived the differential 
equation for H'^{x, 0,p). 

Additional properties and further generalizations where studied by Chat- 
terjea Chongdar [3], Joshi and Prajapat J^, Munot, and Mathur [T7] . 
Rajagopal [iH], Riordan [inj, Saha [20], Shrivastava [22j and Singh and Ti- 
wari |2n]- 

In this work, we will consider the polynomials 

= (-l)"exp (xO :^exp {-x^) , (1) 
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with r = 2,3,..., n = 0,l,..., which correspond to the particular case 
if^(x,0, 1). Clearly H^{x) = Hn{x) = Hermite polynomial of degree n. The 
first few H^{x) are 

H;^{x) = 1, Hl{x) = rx'-\ H;{x) = r^x^^'-^^ - r (r - 1) x^-^ 

Hl{x) = r^x^^'^-'^^ - 3r2 (r - 1) x^'-^ + r (r - 1) (r - 2) x'-^ 

and in general 

= r"x"('^-i) (-r)„ x''-", 

where (■)^ denotes the Pochhammer symbol. Our work was motivated by 
the talk "Asymptotics for Hermite type polynomials", delivered by Professor 
Wolfgang Gawronski at the conference organized in honor of the 65th birth- 
day of Nico M. Temme in Santander, Spain on July 4-6, 2005. In his talk, he 
considered the behavior of H^{x) and its zeros using Plancherel-Rotach type 
asymptotics. 

In Section 2, we present the basic properties of H^{x). Although some 
are not new, we present proofs of all of them for completion purposes. In 
Section 3, we present the asymptotic analysis using a modified ray method, 
developed by Dosdale, Duggan and Morgan in |8j and formalized by Costin 
and Costin in In a previous work [7 , we successfully applied the same 
technique to the classical Hermite polynomials H^{x). Section 4 contains our 
main result and supporting numerical examples. 



2 Properties 



The Hermite polynomials admit the simple hypergeometric representation 



Hr,{x) = (2X)" 2^0 



{n-l) 



X 



where pFq [■] is the hypergeometric function lOj. On the other hand, for 
H^{x) we need to consider the extension of pFg given by Meijer's G-function 
^. 

Proposition 1 The polynomials H^{x) can be represented in terms of Mei- 
jer's G -function by 



H!^{x) = exp {x^ 



x) ^'^'^'+1 



X 



1 - 



0,1 + ^, 1 + ^,' 
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Proof. We have 



exp {—X 



E 

j=0 



dx"^ j\ 



^ oo 

-E 



{r])\ {—x'^y 



x-^{rj-n)\ j\ 



(2) 



Using the multiphcation formula for the Gamma function [21], we obtain 



r-l 



{rjy. r (rj + 1) 



i=0 



(rj-n)! T{rj -n + 1) ^ , 

rn-n+ij]r(^+j) 

i=0 



Replacing in Q we get 



n 



1 V r /V r / 1 



exp [—X 



f \ n 
X 



=1 \ r J j=0 i=l \ r Jj -J' 



r\n^ r (^) 



lip 

4=1 V r / 



1— n 2— ra _ _ _ r—n 



(3) 



Thus, from the definition of H^{x), it follows that 



H^{x) = exp {x''] 



X 



friKiL F 

lip / i-n \ ^ r 
i=l \ r J 



1 2 



,1 



1— n 2—n r—n 



\- r ' r 



X 



(4) 



Using the relation between the hypergeometric function and Meijer's G- 
function 



1=1 



pJ- q 



ai, ■ ■ ■ ttp 



X 



— X 



1 — ai, ■ ■ ■ 1 — ttp 
0,1 -1-6, 



i=l 



in (jH), the results follow. 
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Proposition 2 The polynomials H^{x) satisfy the differential-difference equation 

d 



HnA'^) + Yx^n{x)=rx^-'H:,{x). 



(5) 



Proof. The result follows immediately from the Rodrigues formula since 
d d^ 



d^^ 

dx^' 



+ (-1)" exp -p-^ exp [-X- 



rx^-'H:{x)-H:^,{x). 



When r = 2, we recover the well-known formula for the Hermite polyno- 
mials 

Hn+l{x) + = 2xHn{x). 

Proposition 3 The polynomials H^{x) have the exponential generating function 



G{x, t) = J2^ni^'>ni. ^ • 



n=0 



Proof. From we get 



G'(x,t)=exp(x^)5^(-ir- 



n=0 

oo 



(it" 



exp (x'') V — 

n=0 

exp (x'') exp [— (x — t)*^] . 



J t=0 



(6) 



In particular, for r = 2, we have 

exp [x"^ — (x — t)^] = exp (2tx — t^) , 
which is the exponential generating function of the Hermite polynomials PE] • 



5 



Proposition 4 The polynomials H^{x) admit the explicit representatu 



rk—n 



where 



3=0 



3 J \ n 



Proof. From (HI) we have 



Hni^) = (-1)" exp {x') — exp (-x" 



:-ir 



.fc=0 



A;! 



E 



,i=o 



X" 



oo / „ X fc 

X 



E 

,fc=0 



A;! 



E 

.j=o 



{-ly frj 



[X 



and hence |2ni (2.8)] 



H:(x) 



:-irn\ 



X" 



E 

fc=0 



E 



_iy / I 



j\ \nj{k-j)\ 



[X 



If rk < n, we have 



and therefore, 



0, 0<j<k 



n I CO 



X" 



E 



E 

.j=o 



(-1)^' /r A 1 



\x 



from which (0) follows. 
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Note that when r = 2, ((7j) reduces to the well-known representation of 
the Hermite polynomials jlbj 



Hjx) 



-l)"n! 



Li=o 



/2k-n 



X 



(2A;-n)!(n- A;)! 



2fc-n 



L2J / ,nA; 1 



A:=0 



(n-2A;)!A;! 



Corollary 5 The polynomials H^{x) satisfy 



H:(x)= u'^HLiux), with 



(9) 



The symmetry relation generalizes the reflection formula of the Her- 
mite polynomials [Tl)] 



HJx) 



Remark 6 It follows from (0j that it is enough to analyze H^{x) in the 
region |arg(x)| < ^ of the complex plane. In particular, the roots of H^{x) 
are completely determined by its positive real roots. To show this, in Figure 
U\we plot the roots of II^{x) in the complex plane. 



Proposition 7 The polynomials H'!^{x) can he represented by 



Ni + ---+Nr=nj=l ^-'^ I 

n>0 V 



r\ ^ (-1) J (!— ilAT, 

X i . 



n>0 

In particular, we have 



Kio) = i-iy-^'' j^n^n) , (10) 

where JM (2.32)] 

n,w = iy:exp(H^^U| i„ ;{ --I' . (11) 

^ ^ \ r J 10 otherwise ^ ' 

i=i ^ ^ ^ 
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Figure 1: A plot of the zeros of H^{x) in the complex plane. 
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Proof. From (jU)) we get 



2^- 

n=0 



exp 



exp [x^ — {x — ty] = exp 



X 



j=0 ^"'^ 



JJexp 



j=lAf,=0 



r \ J (—1) J ('■-j) 



^11 



X 3 



A^iH |-A^r=" 

n>0 



i=i 



■1) 



(j-i) 



-X ^ 



Ni 



where we have used equation (2.10) in 

Comparing (fT^ and (fT^ . we conclude that 

Ni+---+Nr=nj=l ' 
n>0 



NiA \-Nr=nj=l 

n>0 



In particular, when x = 0, we must have 

Ni = ■■■ = Nr-i = and A^^ = n, 
from which ()10|) follows. ■ 
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Proposition 8 The polynomials H^{x) satisfy the recurrence relation 

r-l 



K^^{x) = r 



n \ / r — 1 

k 



(14) 



Proof. From (jH)), we have 



dG 
'dt 



r{x-ty-^G 



or [Will], 



n=0 



fc=0 



r — 1 



n=0 



Comparing coefficients of t we get 



n! 



fc=0 ^ 



r — 1 
A; 



.-l-fc-^n-fe(^) 



(n- A;)! 



and (HD) follows. ■ 

When r = 2, (|THl reduces to the three-term recurrence relation for the 
Hermite polynomials |il? 



Hn+iix) = 2xHJx) - 2nHn-^i(x). 



Proposition 9 Let Un{x) = exp {—x^') H]^{x). Then, 



k=0 



fr -l\{n + r - 1)! ^ 



k J {n + k)\ 



-x u. 



(15) 



where 



(k) 



dx^ 



Un (x) . 



Proof. It is clear from (jT)) that 

Un+k{x) = {-if u'^^\ 



(16) 
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From (HD), we get 



Thus, 

/ N v^/ /r - 1\ (n + r - 1)! 

= r > (-1) — ———XUn+k{^)- i^V 



fc=0 ^ ^ 



(n + fc)! 



Using (fTH|) in (fT7|) . the result follows. ■ 
For the case r = 2, (fT3j) gives 

[exp {-x^) Hn]" + 2a; [exp {-x^) + 2 (n + 1) exp (-x^) Hn = 0, 
which is equivalent to the differential equation of the Hermite polynomials 

m 

H'l - 2xH'^ + 2nHn = 0. 



3 Asymptotic analysis 

We seek an approximative solution for (jSj) of the form 

H^{x) exp[f{x,n) + g{x,n)] , n ^ oo (18) 

with 

3 = o{f), n^oo. 
Since Hq{x) = 1, we must have 

/(x,0) = and ^(x,0) = 0. (19) 

Using (fTHj) in (0), we have 

exp + ^ + + ^ + (20) 

\ dn 2 d'n? dn J 

+ ill- + ^ ) exp {f + g) = rx'"^ exp (f + g) , 



dx dx 
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where we have used 

df 1 

/(x, n+1) = f{x, n) + —{x, n) + 2 ^(^' 

From we obtain, to leading order, the eikonal equation 

9f\ , a/ r-l 



and 

^ ^ \ 2 dn J dx ^ \ dn J ' 

or, to leading order, the transport equation 

dg dg f df 



2 drfi ^ dn^ dx^ \ dn' ^ ^ ^ 

3.1 The rays 

To solve ()2H) . we use the method of characteristics, which we briefly review. 
Given the flrst order partial differential equation 

F{x,n,f,p,q) = 0, with p = 9 = 

ox on 

we search for a solution f{x,n) by solving the system of "characteristic 
equations" 

dx dF dn dF 

dt dp ' dt dq ' 

dp _ dF dF dq _ dF dF 

dt dx ^ df dt dn ^ df 

df _ dF dF 

dt ^ dp ^ dq' 

with initial conditions 

F [x(0, s), n(0, s), /(O, s),p(0, s), g(0, s)] = 0, (23) 

and 

^/(O, s) = p(0, 5)^a;(0, s) + g(0, s)-^^n{0, s), (24) 
12 



where we now consider {x, n, f,p, q} to all be functions of the variables t and 
s. 

For the eikonal equation (PT|) . we have 

F (x,n, f,p,q) = e'^ + p — rx'^~^ (25) 
and therefore the characteristic equations are 

dx dn q dp ( r-2 dq 

-dt=^^ -dt='^ -dt='^'-^^'' ' rft=°' 

and 

f^=P + qe'. (27) 
Solving (pUj) subject to the initial conditions 

x{0,s)=s, n{0,s) = 0, q{0,s)=A{s), 

with A{s) to be determined, we obtain 

x = t + s, n = te'^, p = r {t + sY^^ + B{s), q = A, 

for some function B{s). From ()2H|1 we have 

+ rs'-^ + B- rs'-^ = 

and B = — e^. Thus, 

x = t + s, n = te'^, p = r{t + sY~^ -e^, q = A. (28) 
Since (0 implies that /(O, s) = 0, we have from ^ and ^ 
(rs*-^ - e^) X 1 + A X = 0. 
Hence, A = In {rs^~^) and therefore 

x = t + s, n = rs'-H, (29) 

p = r [(t + sy-^ - s'-^] , g = In (rs*-^) , (30) 

with s > 0. 
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Figure 2: A sketch of the rays with r = 5 showing the caustic. 



3.2 The caustic 

Sketching the rays (j^^ . we observe that they fill the region x > Xc{n), 
where Xc{n) is the caustic, i.e., the points in the (x, n)-plane at which the 
Jacobian 

^ , . dx dn dx dn ^ o ^ , , . 

^ a aJ-aJa="K'--i"-'l P^' 

is zero (see Fig|2l). 

From (j!?T|l we have, for s > 0, 



J(t,s) = ^ s = (r - l)t. 
Usine (El in (121 we obtain 



X 



with 



J(t, s)=Q<^t = -<^s = \x 
r 



r - 1 1 

A= , -<A<1. 

r 2 



(32) 

(33) 
(34) 
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From (j221) and (jSSl) we conclude that 

3.3 The functions / and g 
Using dnO)) in dHj) we have 



(35) 



dt 



■ [{t + sy-^ - s'-^] + In {rs'-^) rs'-\ 

while (ITI?|l implies that /(O, s) = 0. Solving we obtain 

/(t, s) = (t + s)" - + [In (rs"-^) - l] rs'-H 

or, using (j^ . 



(x — t)'' + n 



In 



To solve the transport equation (j^ . we need to compute |^ 
as functions of t and s. Use of the chain rule gives 



(36) 

(37) 

(38) 

and |2 

ax 



' dx 


dx ' 




-dt_ 


dt_- 




'1 


0" 


En 






dx 


Qn 






on 




OS 


ds 







1 


.dt 


ds. 




-dx 


dn. 





and hence, 



-dt_ 


dt_- 


1 


dn 


dx' 


dx 


dn 






ds 
-dx 


ds 
dn. 


~ Jit.s) 


dn 
. dt 


dx 
dt . 



where the Jacobian J(t,s) was defined in ^T\i . Using and (jHTj) 
we have 



-§1 dt_ 

dx dn 
ds ds 

. dx dn 

From (jnil) and (gOl) we get 

d\f d^f dt 



[r - l)t 



(r - l)t 
—s 



(39) 

in dnni) 

(40) 



+ 



dv? dndt dn dnds dn rs^' ^[(r — l)t — s]' 



dg dg dt dg ds 

dn dt dn ds dn 

dg dg dt dg ds 

dx dt dx ds dx 



ds 



dg 

dt 



rs^ 2 ^(^j. _ _ 
(r 



(41) 



-l)^lf-^lf 
(r-l)t-s 
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Using (jUD in ((22D, we obtain the ODE 



dg 



r — 1 



dt 2 [(r - l)t - s] ' 
and (fT^ gives g{0,s) = 0. Solving (jl^ we get 



(42) 



g{t,s) = - In 



s - (r - l)t 



or, using (jSH)), 



1 / x - t 

^ = o 7 

2 \x — rt 



(43) 



(44) 



Note that g is undefined when x = rt, i.e., for x = Xc{n) 
Thus, for X > XJn) we have 



X 



(x — ty + nln 



n 



X — t 
X — rt' 



n ^ oo 



(45) 



(46) 



Hl{x) ~ exp 

with t{x, n) defined imphcitly by 

r {x — ty^^ t — n = 0. 

3.4 The function t{x, n) 

To solve PUI) . we shall use Lagrange's inversion formula. 



Theorem 10 Let tl^{u) and (p{u) be formal power series in u, with (p (0) = 1. 
Then there is a unique formal power series u = u{z) that satisfies 



u = Z(f){u). 



Further, we have 



[z'] {^Hz)]} = -[u''-'] {^\u)c/>{uy}, 



(47) 



(48) 



where by [2;'^] {i/j^z)} we mean the coefficient of z^ in the power series of 
ij{z). 
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Proof. See j2Sl Theorem 5.1]. ■ 

We first rearrange ()4(jj) so that it looks like ()47|1 and obtain 



t n 



X rx' 



1 - 



X 



l-r 



Using (jlH|l we then have 



n 

rx^ 

1 

k 



t 

X 

k-l 



k 

oo 

E 

j=0 



(l-r)fe' 



X 



Thus, 



X ^k 



k=l 



1 /(l-r) k 
k[ k-l 



or 



^1 f(l - r) k 
tont{x,n) = x22l[ 
k=i ^ 



where A was defined in flMl) and 



(-1) 



rx^ 



Xx 



1-p 



n 



rx' 



fc=0 ^ ^ 



, (49) 



(50) 



To find the region of the complex plane where p {z) is analytic, we com- 
pute its radius of convergence using the ratio test. We have 



rk 
k 



1 



l—rk 



lim , ^ , 

fc-»oo I r {k + 1) 



r — 1 ' 



where we have used Stirling's formula and 




2^6 ^, Z — > OO. 
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Thus, we conclude that to^t{x,n) is analytic for 



n 



rx"' 



< 



r — 1 



or |x| > Xc{n), with Xc{n) defined in (jHSI)- Note that from we have 

tout{t^x,n) = LL!tont{x,n), for to'' = 1. 
Proposition 11 The function p{z) satisfies the following properties: 
1. We can represent p{z) as a hypergeometric function 



P (z) = r-lFr-2 



^12 

1 2 r-2 



r-2 

r 
r-2 

r-1' r-1' • • • ' r-1 

^. In particular, for r = 2, we have 



jr-l) 
X'- 



z < 



A'" 



r — 1 



(51) 



p{z) = Vl-Az, \z\<- 



and for r = 3, 

p (x, n) = cos 



Proof. 

1. We have 



2 • / 3 /- 

- arcsm -v 

3 V2 



z < 



27' 



-iF, 



r-2 



-112 



7 ' r ' * * * ' 



^ 1 rp 1 rp 

1 2 
r-1 ' r-1 ' 



r-2 



r 
r-2 



r-1 



X- 



oo r-3 



fc = Oj: 



(r-1) 



r-1 

oo 11 \ r Jk f 1\ 



E 

fc=0 



r--2 (r^) (I) 

j=0 



jr-l) 



— ^ r(^ + A;) -^ r(gf) 1 
hl\ r(^) l{ri^^ + k){i-rk)ki 
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k=0 

oo 

E 

k=0 



r (rk + 1) 
r [(r - 1) A; + 1] 



(r-1) 



r-l 



(1 -rfc) A;! 



(r-l) 



rk\ ( X 



k J \r — 1 J (1 — rk) X 



(r-l) 



p(.). 



2. If r = 2, we have 



fc=0 ^ 



2A; 



Using the identity [211 (2.43)] 

f:(2*).' = (i-4xn 



we have 



E 

fc=0 



2A;\ 1 



k J l-2k 



-X 



-2k 



i/E 



2k 



X I ^ — ' \ k 

k=0 



X 



-2k 



I r _l 1 , 

-/ (1-Au-^) ^ du = -^/^^^^^, 

X J ^ ' X 



or 



fc=0 



2k\ 1 



e: T 



2k 



X 



VI -4x 



and ()52|) follows. 

The case r = 3 was computed using (j3T|) and Maple 10. 



(52) 



We shall now find a representation for t(x, n) in the disc |x| < Xcin). We 
first observe that when x = 0, we have from pUj) 

r{-\Y'^\t{^,n)X -n = ^, (53) 

which implies that t(0,n) = r;(n), where 

ri(n) = [n(l-A)]^"^exp[A(2/ + l)7ri], 1 < / < r. (54) 
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Using in (jiUj) . we get 



1, 



or 

To solve we use the following Lemma. 
Lemma 12 Given an algebraic equation of the form 



\ r-l 

= 1. (55) 



a{a-bY'^ = 1, Cy^O, (56) 



we formally have 

fc=0 ^ 



Proof. Solving for 6 in (|56p we get 

1 

b = a — a^-" . 

Letting = a^^-i — 1, we have a = (^ + 1) and therefore, 

6 = (e+i)'^-(e+i)"' = e(e + i)"' 

or 

e = 6(e+i)'. (57) 

Applying Theorem ITUl to (IKTj) with (^ = (^ + 1)^ and ^ (^ = (^ + 1)^ , 
we obtain 

[fc'^] {a (6)} = [b'] [e (6)]} = I [e-'] (e + 1)"' (e + 1)' } 

and the result follows. ■ 
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Thus, applying the Lemma to (fH3|) we find that 



j=0 



j \ j J \n) 



which we can write as 



tin(x, n; I) = Ti{n) + Az/x 



X 



with 



Ti{n) 
1 



1 < / < r, 



k k+l 



Applying the ratio test to (jSHl), we get 



lim ■ 



k J 



fc+i \ 

r \ 1 

A; + 1 



A-^ (1 - A) 



A-l 



il-A 



(58) 
(59) 



From ()54p we have 

\r,{n)\ = [n{l-\)y 
and therefore tiri{x^n] I) is analytic in the sector 

|x| < [n (1 - A)]'-^ A-^ (1 - A)^-' = n^-^A-^ = X,{n). 

The function /i {z) doesn't have a simple expression in terms of elementary 
functions, except for r = 2, when we have 



X 



2 + V4 + 2;2 



4 Summary and numerical results 

We summarize the results of the previous section in the following theorem. 
Although our analysis was done on the positive real axis, we can extend our 
results to the whole complex plane, with the exception of the caustic circle 
\x\ = XJn) = X~^n^~^. 
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Theorem 13 Let 



Kix) = (-1)" exp (x*^) — exp (-x' 



Then, as n oo, 



1. For \x\ > Xc{n), we have 



H'^{x) ~ Hout{x, n) = exp (x' \ 1 



1 - A + Ap 



n 



rx' 



— n 



n 



X 



i-Mi-p(^) 



A = ^— ^ and 



rk 



2. For \x\ < Xc{n), the function t{x,n) is multivalued and therefore we 
need to add all the different contributions, 



H'!^{x) ~ Hin{x, n) = exp ix^—<x — ri{n) — Xxfi 



X 












Ji{n)_ 







X 



n 



Ti{n) + Xx^i 



Ti(n) 



n{n) - 


X ' 




X 

_n{n) _ 




\ 


rTi{n) — X ' 




- (r - 


X 




} 



with 
and 



Ti{n) = - A)]^"^exp [A (2/ + 1) vri] 



fc=0 

In Figure El we sketch the ratio i^^f^^^ , for r = 5, n = 5 and Xc{5) ~ 
1.649 < \x\ < 7. We clearly see how the approximation breaks down inside 
the caustic region |a;| < Xc{n). 

In Figure Ufa) we compare the values of H^{x) and Hi^{x,n) for r = 5, 
n = 5 and |a;| < Xc{5). To show in detail the graphs for values of x close to 
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Figure 3: A sketch of ^ ^(5) the outer region Xc{5) < \x 
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Figure 5: A comparison of H^{x) (solid curve) and Hont{x,5) (ooo) with 
X = 2e'^. 

the zeros of H^{x), we sketch the two functions in smaller intervals in Figures 
Hb) Mc) andUd). 

As we mention at the beginning of this section, our approximations are 
valid in the complex plane minus the circle = Xc{n). To illustrate this, 
in Figure El we graph the real (a) and imaginary (b) parts of H^{x) and 
Hout{x, n) for r = 5, n = 5 and x = 2e'^. From Corollary El we know that it is 
sufficient to consider the sector \9\ < |. Finally, in Figure El we do the same 
for the functions H^{x) and Hont{x, n), with r = 5, n = 5 and x = \e^^ . 
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